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By 
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ABSTRACT 

This is the second of six volumes, each bearing the same report number, 
but dealing with separate problem areas concerning the stability of 
eccentrically stiffened circular cylinders. The complete set of documents 
was prepared under NASA Contract NAS8-11181. This particular volume presents 
a practical method for the determination of design values for the critical 
buckling stresses of axially compressed, curved isotropic skin panels having 
various degrees of edge restraint. For wide panels which behave essentially 
as full cylinders, the method provides lower-bound predictions. As the 
panel width decreases, a transitional relationship is employed which, for 
narrow panels, approaches standard flat plate values. In order to facilitate 
ready application, this volume includes several groups of design curves as 
well as a digital computer program which provides either automatic plotting 
or single-point solutions, as desired. 
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DEFINITION OF SYMBOLS 


Symbol 


Definition 


a Axial length of curved isotropic skin panel; 

Spacing between rings. 

b Circumferential arc length of curved isotropic 

skin panel; Spacing between stringers. 

C Factor defined by equation (2-8). 

E Young's modulus. 

K Buckling coefficient for a flat plate. 

K Buckling coefficient for an isotropic flat plate 

c having the loaded edges simply supported and the 

unloaded edges fully clamped. 

K Buckling coefficient for an isotropic flat plate 

s having all edges simply supported. 

R Radius to middle surface of basic cylindrical skin. 

t Skin thickness. 


v Poisson's ratio. 

0 CL Classical theoretical value for the compressive 

buckling stress of an isotropic cylinder ( = .606 -g- ). 

a Critical compressive stress, 

cr 

a Critical compressive stress for the buckling of 

P an isotropic flat plate. 

o Critical compressive stress for the buckling of 

an isotropic cylinder. 
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SECTION 1 
INTRODUCTION 

This is the second of six volumes, each bearing the same report 
number, but dealing with separate problem areas concerning the stability 
of eccentrically stiffened circular cylinders. The complete set of 
volumes was prepared under NASA Contract NAS8-11161. This particular 
volume deals with the buckling of the curved isotropic skin panels which 
are bounded by stringers and rings. Although the buckling of these skin 
panels is usually not catastrophic, the related critical stress values must 
be established in order to compute effective skin widths for use in the 
investigation of other modes of instability such as so-called panel in- 
stability (see GLOSSARY, Volume I [l])and general instability (see GLOSSARY, 
Volume I Cll) . In addition, although the overall stiffened cylinder can 
usually continue to support loading beyond the point at which the isotropic 
skins buckle, this condition will sometimes be undesirable from other view- 
points (fatigue, excessive deformation, etc.) 


Numbers in brackets [ 1 in the text denote references 


listed in SECTION 6. 
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SECTION 2 
EQUATIONS 

2 . 1 GENERAL 

The procedures of this volume are based on the Schapitz [ 2 ] criterion 
for the buckling of curved isotropic skin panels. This criterion supplies 
a means for evaluating the effects of skin panel geometry as the trans- 
ition is made from wide panels behaving essentially as full cylinders to 
smaller panels which approach the behavior of flat plates. In the former 
case the fixity effects at the longitudinal edges are negligible. However, 
for the narrow skin panels this fixity has a marked influence on the critical 
buckling stress. 

In the application of the Schapitz criterion, three different empirical 
foundations were considered for the case where the longitudinal stiffeners 
are widely spaced. One of these results in typical strength predictions 
which would primarily be of use in making comparisons against test data. 

The other two give design levels of high reliability recognizing the wide 
scatter which exists in actual strengths. The recommended design procedure 
was selected from these latter two alternatives and is consistent with the 
isotropic cylinder criterion of reference 3. 

It is helpful to note that the character of flat plate buckling is 
much different from that exhibited by wide curved skin panels and isotropic 
cylinders. The flat plate can continue to support steadily increasing 
external loads well into the postbuckling region. This is in contrast to 
the sudden drop-off in load usually observed for wide skin panels and 
isotropic cylinders. Consequently, the Schapitz criterion [ 2 ] utilizes 
full theoretical predictions as the boundary case of a flat plate is 
approached. Therefore, within this region, one might expect that the test 
data will display some small degree of scatter on both sides of the 
proposed design curves. However, because of the physical behavior cited 
above, this generally will not lead to any serious structural deficiencies. 
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2.2 SCHAPITZ CRITERION 

The Schapitz criterion [ 2 ] is embodied in the following simple 
relationships which have been verified by the rederivation of reference 4: 

( 2 - 1 ) 
(2-2) 

(2-3) 

then a = a (2-4) 

cr K 

where , 

o = Critical stress for buckling of a flat 
p isotropic skin panel. 

c R = Critical stress for buckling of an 

isotropic cylindrical shell. 

o = Critical stress for buckling of a curved 

cr isotropic skin panel. 

Equation (2-2) supplies the transitional relationship. To clarify dis- 
cussions which follow, it is pointed out here that the notation of Figure 1 
will be employed in the application of the foregoing criterion. Note that 
the dimension b is the true circumferential arc length of the curved panel. 
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2.2.1 FLAT PLATE EQUATION 


From equations (2-1 ) through (2-4), it can be seen that application 
of the Schapitz criterion t2] requires the availability of suitable means 
for computation of the flat plate critical buckling stress . The 
following well-known formula [5] will be used for this purpose: 


cr 


P 


K 


E 

(1-v 2 ) 



(2-5) 


where 

K = Flat plate buckling coefficient which 
depends on aspect ratio (a/b), boundary 
conditions, and type of loading. 

E = Young's modulus, 
v = Poisson's ratio. 


In this volume it will always be conservatively assumed that the loaded 
edges are simply supported. No advantage will be taken of any rotational 
restraint along these particular boundaries. However, for the longitudinal 
edges, both the simply supported case and the fully clamped condition will 
be covered. Therefore, the following notation may be used to identify the 
buckling coefficients K of interest here: 

K = Buckling coefficient for flat plate having 
s 

all four edges simply supported. 

K c = Buckling coefficient for flat plate having 
the loaded edges simply supported and the 
remaining two edges fully clamped. 

Thus, for the purposes of this volume, it is convenient to rewrite equation 
(2-5) in the form 

a = K r -I ^o-(r) 2 (2-6) 

p sLc:i d-v 2 )V b ' 

to identify the nature of the selected boundary conditions. 
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2.2.2 FULL CYLINDER CRITERIA 


From equations (2—1 ) through (2-4), it can be seen that application 
of the Schapitz criterion [2] also requires the availability of suitable 
means for computation of the critical buckling stress of an isotropic 
cylinder. It must be kept in mind, however, that the test data for 
axially compressed isotropic cylinders usually fall far below classical 
theoretical values and display a high degree of scatter. Therefore, design 
stress levels for such cylinders are normally established by applying 
empirical reductions to the classical theoretical predictions. In this 
volume, three different empirical criteria are considered. These al- 
ternatives are identified below as OPTIONS 1, 2, and 3. 


OPTION 1) 


The lower-bound approach of Seide, et al. C 6 1 : 



where 


C 


= 0.605 - 0.546 




(2-7) 

( 2 - 8 ) 


OPTION 2) The ’'best-fit'* (or "mean-expected") relations of reference 7 
which may be rewritten in terms of the parameters of interest 
here as follows: 


— Jl < i o 
b R A * 


SI* ^ 


“-“(I) 1 ' 6 + °- 109E (f SI)' 

C R = 11.28 (f)' 1 ' 6 *0.109 e(S f |J 

- 1 - 418 (f)' 1 ' 6 ***. ( I I) 


(2-9) 


( 2 - 10 ) 
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OPTION 3) The statistical level C 7] which represents 90% 
probability with 95% confidence; i.e., there is 
95% confidence that 90% of specimens tested would 
exceed the buckling levels given by: 

ll <1 - 0 *° R = 8.° lE (f)" 1 ' 6 * 0.076E (f {t |)’ 1,3 (2-11) 

! 1 * I- 0 - 7 - 52E (f) " 1 ' # * °- 07 “ (f f |) * 1,3 

-!.418 (I)' 1 ’ 6 Lo g> (2-12) 

In view of its relative simplicity and apparent reliability (see 
SECTION 3, "TEST DATA COMPARISONS"), it is recommended here that OPTION 1 
be used in the analysis of curved isotropic skin panels bounded by stringers 
and rings. However, it is further recommended that the OPTION 1 application 
be restricted to cases where (a/R) s 1. Since practical configurations will 
almost always satisfy this condition, it does not constitute a severe 
limitation. This restriction has been imposed here in view of the situation 
depicted in Figure 2 where the lower-bound curve of Seide, et al. C 6 3 is 
compared against a 90% probability curve (a/R = 1.0) of reference 7. A 
comparison could be made against the corresponding 99% probability curve of 
reference 7 but it is thought here that this latter statistical level is 
excessively conservative for multiple load-path structures such as the 
stiffened cylinders under consideration. The curves for 99% probability are 
more suitable to single load-path configurations such as pure monocoques. 

The comparison given in Figure 2 shows that the curve of Seide, et al. [6] 
corresponds reasonably well to the reference 7 curve for (a/R) =1. The 
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1.0 


cr 


'CL 



(?) 

Figure 2 - Comparison of Lower-Bound 
Curve of Seide, et al. [6l 
Against 90% Probability 
Curve of Reference 7 
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Seide criterion is somewhat more conservative in the lower (R/t) region. 

On the other hand, for (R/t) values in excess of 1000, the Seide criterion 

is the less conservative of the two alternatives but in this high (R/t) region 

, 2 

the differences in terms of lbs/in. are of rather small magnitudes. 

It should be noted that, unlike OPTION 1, the equations for OPTIONS 2 
and 3 include a statistical effect from the cylinder length. This added 
influence complicates the graphical presentation of the Schapitz criterion 
[2] to such a degree that the associated buckling curves are considered here 
only as supplements to the recommended curves given in SECTION 4, "DESIGN 
CURVES". These supplementary plots are furnished in Appendix A. The 
buckling curves based on OPTION 2 could be of some use in the evaluation 
of test results while the curves based on OPTION 3 might be used to obtain 
design values for cases where (a/R) > 1. 

2.2.3 GRAPHICAL REPRESENTATION 

In order to permit a graphical representation of the buckling criterion 
of equations (2-1) through (2-4), equation (2-6) may be rewritten as 


a 

P 



(1 -'' 2) (ff(l ) 2 


(2-13) 


Within a nondimensional, logarithmic format, equations (2-1) through (2-4) 

can then be presented as shown in Figure 3. The transition curve defined 

by equation (2-2) becomes tangent to the full cylinder curve when 

©„ = 2o . For all (R/t) values greater than that of the tangency point, 

R p 

the skin panel behaves as a complete isotropic cylinder. For all other 
(R/t) values, the transitional relationship applies. Note that the trans- 
ition curve asymptotically approaches the curve for c . 

r 
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H 


Figure 5 - Nondime nsional Logarithmic Plot of 
Buckling Criteria for Isotropic Skin Panels 
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SECTION 3 

TEST DATA COMPARISONS 

In order to evaluate the reliability of the recommended methods 
of this volume , comparisons were made between predicted critical stresses 
and experimental data appearing in the literature. In particular, the 
test data of references 8, 9, 10, and 11 were compared against Schapitz 
criterion [2] predictions which employ the OPTION 1 full cylinder equation. 

The results obtained from this investigation are summarized in Table I. 

All predicted values were obtained using digital computer program 3875 
which is given in SECTION 5 of this volume. It should be recalled that 
the OPTION 1 equation describes the lower-bound curve of Seide, et al. [6] 
for 0^ [see equations (2-1 ) through (2-4)]. Hence one would expect all 
the predicted critical stresses to be conservative with respect to the 
actual test results, particularly in the transition and full cylinder 
regions of behavior. Considerable scatter would also be expected of the 
test data in these regions so that the degree of conservatism would vary 
rather widely from specimen to specimen. Inspection of Table I shows that, 
for most of the data considered, the predictions (based on OPTION l) were 
indeed conservative and that considerable scatter did exist in the test 
results. However out of a total of 127 test points, 24 fell below the 
supposedly conservative predictions. On the surface, this would appear 
to be an unsatisfactory condition and would tend to discredit the re- 
liability of the proposed analysis methods. However, further scrutiny of 
the test data shows that, except for 4 points, the troublesome data came 
from specimens which fell into the flat plate region of behavior. There- 
fore, in general, the discrepant prediction values could only be significantly 
lowered by applying a reduction factor to the classical flat plate values. 

This would be contrary to time-tested current design and analysis practices 
which have usually proven to be perfectly adequate. It has therefore been 
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TABLE X - Buckling of Isotropic Skin Panels 
(Continued) 
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simply describes a lower-bound curve to the actual test data. Only the second-series data of reference 11 were 
used here. The first-series data of that reference are of questionahle validity due to test procedure difficulties. 
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concluded that, in general, the discrepancies under discussion here 
are probably due to the inability to select precisely correct values 
for the buckling coefficient K [see equation (2-5)]. This coefficient 
is strongly dependent upon the actual degree of restraint afforded to 
the skin panel by the boundary members. Further uncertainty is introduced 
through the usual experimental difficulties surrounding the decision as to 
when initial buckling has in fact occurred in an imperfect flat plate. 

For the test data of reference 8, the K value was taken equal to 
5.1. This constitutes an average value obtained from the curves developed 
by Peterson and Whitley [ 12] for long skin panels having Z-shaped stringers 
as longitudinal boundary members. In studying the curved specimens of 
reference 9, use was made of the flat plate data which was concurrently 
obtained for control-type purposes. These flat plates were subjected to 
boundary conditions which were similar to those of the corresponding curved 
skin panels. The experimental critical stress for each flat plate was 
substituted into the equation 


K = 


a 

P 


(1-v 2 ) 

E 



(3-1) 


which is obtained from a simple transformation of equation (2-5). The 
average of these several values was computed to be 5.4. As noted in 
Table I, this furnished the basis for analysis of the reference 9 curved 
skin panels. This same procedure was followed for the specimens of reference 
10. In this case, the average K value was found to be 4.4. As stated 
in the footnote to Table I, reference 11 did not include any listing of 
specific test results for the individual specimens. However, sufficient 
information was provided to ascertain that an average K value of 6.0 
would be applicable there. Hence the related predictions of Table I were 
based upon this value. 

In order to further clarify the preceding discussion, the qualitative 
presentation of Figure 4 is offered. This figure shows the general 
characteristics and relative positioning for each of the following when 
transformed into the desired nondimensional logarithmic format: 
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Figure 4 - Schematic Nondimens ional Logarithmic 
Plot of Test Data for Carved 
Isotropic Skin Panels 
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(a) The flat plate relationship given as equation (2-5). 

(b) The recommended design curve based upon the substitution 

of OPTION 1 cylinder expressions into the Schapitz criterion 

[ 2 ]. 

(c) The classical small-deflection buckling equation for isotropic 
cylinders. 

Also shown in Figure 4 are the individual approximate locations of all 
test points listed in Table I. In the course of the study under NASA 
Contract NAS8-11181, quantitative plots were made for each specimen and 
the corresponding test point was accurately located on the appropriate 
plot. Based on these many different plots, the test points were inserted 
in Figure 4 in approximation to their actual position relative to the 
several basic curves and regions of behavior. In general, this figure 
shows that 

(a) Except for the region of flat plate behavior, the test data 

fall above the recommended design curve but below the value 
which would have been predicted if were computed from the 

classical small-deflection equation, 

(b) All but 4 of the unconservative predictions occur in the region 
of flat plate behavior. 

In view of the foregoing^ it is concluded that the use of OPTION 1 equations 
in conjunction with the Schapitz criterion provides reliable critical stress 
values which can be safely employed for the design of isotropic skin panels 
in stiffened circular cylinders. 
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SECTION 4 
DESIGN CURVES 

All of the buckling curves presented in Figures 5(a) through 5(1) 
are based upon the use of the OPTION 1 full cylinder equations [see 
equations (2-7) and (2-8)] in the Schapitz criterion [ 2 ]. These curves 
were obtained by using digital computer program numbered 3875 and an 
automatic plotter. This program is described in SECTION 5 and can be used 
to obtain additional plots and/or single-point solutions, as desired. 

Since the automatic plotting machine does not provide any capability for 
the print-out of lower case letters, the ratio (a/b) appears on all the 
given buckling curves as (A/B). 

It should be recalled that the Schapitz criterion [ 2 ] considers two 
regions of response dependent upon the two values and o^. When s 2c^, 
the critical stress for the curved skin panel is established by a relation- 
ship which provides a smooth transition between flat plate theory and the 
full cylinder behavior. For > 2c^ the critical stress is simply taken 
as a cr = a^. To establish the point of separation between these two regions, 
this section includes curves 6(a) through 6(d) which show b/R versus R/t 
for the condition o„ = 2 o . Here again, the plotting machine converted 

“ p 

lower case notation into upper case print-outs. All of the curves in 
Figure 6 are based on the OPTION 1 full cylinder criterion. 

The buckling curves of this section are to be applied only where the 
condition (a/R) S 1 is satisfied. For longer skin panels, unconservative 
predictions can result. Whenever (a/R) > 1, it is recommended that the 
predicted critical stresses be based upon the use of the OPTION 3 full 
cylinder equation in conjunction with the Schapitz criterion [ 2 ], A 
selection of design curves based on this option are given in the appendix. 

All of the design curves given in this volume apply only where the 
behavior is elastic. Conventional plasticity reduction factors should be 
used whenever the stresses exceed the elastic limit. In addition, all 
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curves are based on a Poisson's ratio of 0.30. The curves for aluminum 
employ E = 10x10° psi while the curves for steel use E = 30x10 psi. 

In order to use the design curves of this volume, one must first 
select an appropriate value for the buckling coefficient K. Table II 
lists such values for the two boundary conditions of interest here. 
Tables III and IV list the various families of curves provided in this 
section. 


TABLE II - Buckling Coefficients For 
Isotropic Skin Panels Subjected 
To Edge Compression 


BOUNDARY CONDITIONS 

(bj 

K 

Unloaded 

Edges 

Loaded 

Edges 

Clamped 

Simply Supported 

0.4 

K = 7.0 
c 

Clamped 

Simply Supported 

0.6 -» » 

K = 5.7 
c 

Simply Supported 

Simply Supported 

0.4 

K = 7.0 
s 

Simply Supported 

Simply Supported 

0.6 

K . 4.0 

s 

Simply Supported 

Simply Supported 

0.8 - » 

K = 3*29 
s 
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TABLE III - Table of Contents for the 
Design Curves "Buckling of 
Isotropic Panels” 


Figure 

Number 

Ordinate 

Abscissa 

(a/b) 


_K_ 

Page 

5(a) 

/ Buckling Stress \ 


0.4 


7,0 

4-4 

\Elastic Modulus / 

W 




5(b) 

»» 

II 

0.6 - 

00 

5.7 

4-5 

5(c) 

it 

tt 

0.6 


4.0 

4-6 

5(d) 

ii 

tt 

0.8 - 

00 

3.29 

4-7 

5(e) 

Buckling Stress-psi 
For Aluminum 

(+) 

0.4 


7.0 

4-8 

5(f) 

tt 

II 

0.6 - 

00 

5.7 

4-9 

5(g) 

n 

II 

0.6 


4.0 

4-10 

5(h) 

it 

It 

0.8 - 

00 

3.29 

4-11 

5(i) 

Buckling Stress-psi 
For Steel 

(-f) 

0.4 


7.0 

4-12 

5(j) 

ii 

v •• ' 

0.6 - 

00 

5.7 

4-13 

5(k) 

it 

n 

0.6 


4.0 

4-14 

5(1) 

•t 

•i 

0.8 - 

CO 

3.29 

4-15 




TABLE IV 

- Table of Contents for 

the 





Design Curves 
R/t for o R = 

M b/R vs. 
2a " 

P 



Figure 

Number 


Ordinate 

Abscissa 

(a/b) 

K 

Page 

6(a) 

b 

R 

: (For = 20 ) 

R p 

(?) 

0.4 

7.0 

4-16 

6(b) 


It 

tt 

0.6 - 00 

5.7 

4-17 

6(c) 


II 

II 

0.6 

4.0 

4-18 

6(d) 


It 

II 

0.8 -* « 

3.29 

4-19 
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BUCKLING OF ISOTROPIC PANELS 
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Figure 5(e) 
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BUCKLING OF ISOTROPIC PANELS 
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100000 
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Figure 5(i) 
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BUCKLING OF ISOTROPIC PANELS 
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Figure 5(.i) 
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SECTION 5 

DIGITAL COMPUTER PROGRAM 

This section presents the essential features of General Dynamics 
Convair digital computer program numbered 3875. This program was de- 
veloped for the analysis of compressive buckling in curved isotropic 
skin panels which are bounded by the stringers and rings of stiffened 
cylinders. The solution is based upon the theoretical and empirical 
considerations presented in SECTION 2. The reliability of the output 
values is demonstrated by the test data comparisons given in SECTION 3. 

The output can be obtained in the form of automatically plotted buckling 
curves or single-point solutions, as desired. All of the buckling curves 
presented in this volume were obtained by using the automatic plotting 
option of this program. In addition to its primary application to 
curved plates, the program also includes Capability for the analysis of 
flat plates and/or complete isotropic cylindrical shells. 

The input format is shown in Figure 7 . Symbols are listed in Table V . 

A detailed, card-by-card discussion of the input follows below. Runs 
may be stacked. 

CARD TYPE 1: One card per run. Each entry must be integer (s) right 

adjusted in appropriate columns. 

Enter NAB (the number of a/b ratios to be entered on 
cards TYPE 2) in columns 1-5 (30 maximum). 

Enter NBR (the number of b/R ratios to be entered on 
cards TYPE 3) in columns 6-10 (20 maximum) . 
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Enter NBT (the number of b/t ratios to be entered 
on cards TYPE 6) in columns 11-15 (20 maximum). 
Enter NRTIN (the number of R/t ratios to be entered 
on cards TYPE 4) in columns 16-20 (99 maximum). 
Enter NOP (determines which OPTIONS are used) in 
column 25: 


OPTIONS 

NOP 

1 

4 

2 

2 

3 

1 

1 and 2 

6 

1 and 3 

5 

2 and 3 

3 

l,2,an4 „,..3 

7 


As discussed in Section 2.2.2, OPTION 1 utilizes 
the lower bound relation of Seide, et al. [6] for 
d ^ and would be generally employed for design. OPTION 
2 uses mean expected results for and OPTION 3 uses 
90% probability, 95% confidence values for . If 
more than one OPTION is called for (NOP = 3,5,6, or 7), 
complete sets of output for each desired OPTION are 
obtained. 
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Enter IBUG (debugging and printout option) in column 
30. If IBUG as 0 (or is blank), only calculations for 
the R/t values input on cards TYPE 2 will be printed 
out. If IBUG as 1 (or 7*0), calculations will be run 
and printed for 101 values of R/t generated within the 
program whose logarithms are equally spaced. 

Enter K1 (plotting option) in column 35. If Kl = 0 
(or is blank), calculations will be run and logarithmic 
plots will be made using 101 values of R/t generated 
within the program whose logarithms are equally spaced* 
If Kl = 1, no plots will be made. 

Enter K2 (header card option) in columns 36-40. This 
option allows the use of one TYPE 1 header card if 
applicable to a number of runs. If K2 is blank, 0, or 
a TYPE 1 header card will be used as usual. If K2 > 1, 
say "n", then "n" runs may be stacked without repeating 
the header card TYPE 1. When "n" runs have been read, 
the program will then begin reading TYPE 1 cards, if 
additional runs are present. 

Disregard K3 (not used) • 

CARD TYPE 2: There will be NAB cards per run. 

"til 

Enter AB^ (i value of the panel length/width ratio 
a/b to be input) in columns 1-10 (E10.5). 
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CARD TYPE 3: 


CARD TYPE 4: 


CARD TYPE 5: 


thi 

Enter CKNM^ (heading for i output listing and/or 
i*k plot) right adjusted ending in column 20: 

Hollerith characters KS= or KC= or KSsKC= (maximum 
of 6 characters) . 

Enter CAPK. (i tb value of KS or KC) in columns 21-30 

l 

(E10.5). Note that these three entries must correspond; 

x i_ 

i.e., the i TYPE 2 card shows a/b and K or K values 

S C 

which are related since K and K are functions of a/b 
[see equation (2-13)]* 

There will be NBR/6 (rounded to higher whole number) 
cards per run. 

Enter BR values (panel width to radius ratio b/R, 6 to 
a card (6E10.5). 

There will be NRTIN/6 (rounded to higher whole number) 
cards per run. 

Enter RT values (panel radius to thickness ratio R/t), 

6 to a card (6E10.5). 

One card per run. (If K1 = 1 omit this card) 

Enter YBSIG (lower limit of O cr for plots of Buckling 
Stress, psi vs. R/t) in columns 1-10 (E10.5). 

Enter YTSIG (upper limit of O for plots of Buckling 

Ci 

Stress, psi vs. R/t) in columns 11-20 (E10.5). 

Disregard YBBR (not used). 

Disregard YTBR (not used). 

Enter YBSIGE (lower limit of ° cr /E for plots of Buckling 
Stress/E vs. R/t) in columns 41-50 (E10.5). 
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Enter YTSIGE (upper limit of o /E for plots of Buckling 

cr 

Stress/E vs. R/t) in columns 51-60 (E1Q.5). 

CARO TYPE 6: There will be NBT/6 (rounded to higher whole number) 

cards per run. 

Enter BT values (panel width to thickness ratio b/t for 
flat panels only), 6 to a card (6E10.5). Omit this card 
if NBT = 0. These values will be used only if b/R = 0 
( flat plate) . 

A sample input coding form is shown in Figure 8. 

Approximate execution times at the General Dynamics Convair 7094 
DCS installation were 1.5 seconds per a/b and b/R combination for 
NOP = 7 and Kl = 0 . 

The program output consists of a listing of all input in addition 

to calculated results. The calculated results (o /E, o for E = 

cr cr 

6 6 

10x10 psi and o for E = 30x10 psi) are tabulated (and/or plotted) 
for the input values of R/t and the 101 internally determined values of 
R/t, if desired. Separate tables (and/or plots) are formed for each 
a/b and K s £ c «j , b/R, and OPTION combination selected on input. A 
column having the heading "BASIS" is also printed. The symbol I in 
this column indicates that the buckling stress was determined by 
equation (2-2) of Section 2.2. The symbol II in this column indicates 
that the buckling stress was determined by equation (2-4) of Section 2.2. 
The latter case indicates behavior of the panel to be similar to that of a 
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complete cylinder of the same radius, thickness, and length. 

A sample output listing is given in Figure 9 where the first and 
two typical pages are shown for the sample problem of Figure 8. 

Output plots from this sample solution- are included in SECTION 4. A 
basic flow diagram for the program is given in Figure 10. Table VI 
presents a Fortran listing of the program, including input data from the 
sample problem. 
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TABLE V - Program 3875 Notation 


PROGRAM 

NOTATION 

REPORT 

NOTATION 


DESCRIPTION 

AB 

a/b 


Panel length/width ratio. 

NAB 

- 


Number of a/b values input. 

BR 

b/R 


Panel width/radius ratio. 

NBR 

- 


Number of b/R values input. 

BT 

b/t 


Panel width/thickness ratio. 

NBT 

- 


Number of b/t values input. 

CAPK 

K or K 
s c 


Buckling coefficient. 

El , E2 

E 


Elastic modulus, 10x10** and 
30x10° psi respectively. 

FMU 

V 


Poisson's Ratio, v = 0.3 . 

NOP 

- 


Option No. (1, 2, or 3; or 4 
meaning 1, 2, and 3). 

RT 

R/t 


Radius/thickness ratio. 

NRTIN 

- 


Number of R/t values input. 

NRT 

- 


Total number of R/t values, 
both input and calculated. 

SIGRE 

V E 


Nondimensional buckling parameter. 

SIGPE 

V E 


Nondimens ional buckling parameter 

SIGCRE 

a /E 
cr 


Nondimensional buckling parameter 

SIGCR1 

o (for 

Cr R 

E = 10x10 

psi) 

Buckling stress (psi) for aluminum 

fa 

where E = 10x10° psi. 

SIGCR2 

0 (for 

cr 

E = 30x10® 

psi) 

Buckling stress (psi) for steel 

fa 

where E = 30x10 psi. 
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PROGRAM 

NOTATION 

YBSIGE 

YTSIGE 

YBSIG 

YTSIG 

IB 


TABLE V - Program 5875 Notation (Coat'd.) 


REPORT 

NOTATION DESCRIPTION 


Lower limit of ° cr /E for plots. 

Upper limit of a cr / E for plots. 

Lower limit of SIGCR1 and SIGCR2 for 
plots ,psi. 

Upper limit of SIGCR1 and SIGCR2 for 
plots, psi. 

Basis number (1 or 2); e.g., 

BASIS(IB) = BASIS II for 

IB = 2, and c = o . 

’ cr R 
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FORTRAN CODING AND DATA FORM 
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Figure 8 - Sample Input Data-Program 3875 
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LOCAL BUCKLING OF ISOTROPIC PANELS 


NAB = A NBR = 7 NBT = 0 NRTIN 


2 NOP = 4 I BUG 


K1 = -0 K2 


-0 K3 = -0 


PANEL ASPECT RATIO = 0.400 
PANEL ASPECT RATIO = 0.6G0 
PANEL ASPECT RATIO = 0.600 
PANEL ASPECT RATIO = 0.800 


KS=KC= 7.000 

KS= 4.000 

KC= 5.700 

KS= 3.290 


O.CICOOO 

0.020000 

PANEL WIDTH / RADIUS 
0.050000 

0.100000 

0.200000 

0.500000 

1.000000 





RADIUS / THICKNESS 


ICO. 00 


100C.C0 


Y8SIG= 1 • COE 02 YTSIG= 


l.OOE 05 


YBBR= -0. YTBR= -0. 


YBSIGE= 


1.00E-05 YTSIGE= 


I.00E-02 


Figure 9 ~ Sample Output Li sting- Program 3875 
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OPTION 1 


PANEL ASPECT RATIO= 0.80 


KS= 3.290 


PANEL UilDTH/RADIUS= 0.2000 


RACIUS/THiCKNESS 


BUCKLING STKESS/E 


BUCKLING STRESS, PSI 

E = 10,000,000 


BUCKLING STRESS, PSI BASIS 

E = 30,000,000 



IOC. 


9.3 797E-03 

9.3797E 

04 

2.8139E 

05 

I 


1CC0. 


1.4C54E-04 

1.4054E 

03 

4.2161E 

03 

I 


100. 


9.3797E— 03 

9.3797E 

04 

2.8139E 

05 

I 


105. 


8.5763E-03 

8.5763E 

04 

2.5729E 

05 

I 


lie. 


7.8428E-03 

7.8428E 

04 

2.3529E 

05 

I 


113. 


7.1 733E-03 

7.1733E 

04 

2.1520E 

05 

I 


12C • 


6.5619E-03 

6.5619E 

04 

1.9686E 

05 

I 


126. 


6.0037E-03 

6.C037E 

04 

1.8011E 

05 

I 


132. 


5.4940E-03 

5.4940E 

04 

1.6482E 

05 

I 


136. 


5.0284E-03 

5.0284E 

04 

1.5085E 

05 

I 


145. 


4.6032E-03 

4.6032E 

04 

1.3810E 

05 

I 


151. 


4.2148E-C3 

4.2148E 

04 

1.2644E 

05 

I 


15a. 


3.8600E-C3 

3.86C0E 

04 

1.1580E 

05 

I 


166. 


3.5357E-C3 

3.5357E 

04 

1.0607E 

05 

I 


174. 


3.2394E-C3 

3.2394E 

04 

9.7183E 

04 

I 


182. 


2.9686E-C3 

2.9686E 

C4 

8.9059E 

04 

I 


191. 


2.721lfc-03 

2 • 72 1 IE 

04 

8. 1632E 

04 

I 


200 . 


2.4948C-03 

2.4948E 

04 

7.4843E 

04 

I 


209. 


2.2878E-C3 

2.2878E 

04 

6.8634E 

04 

I 


219. 


2.0985E— 03 

2.0985E 

04 

6.2956E 

04 

I 


229. 


1.9254E-C3 

1.9254E 

09 

5.7762E 

04 

I 


240. 


1.7670C-03 

1.7670E 

04 

5.3010E 

04 

l 


251. 


1.6220E-03 

1.6220E 

04 

4.8661E 

04 

I 

01 

263. 


1.4893E-03 

1.4893E 

04 

4.4680E 

04 

I 

1 

275. 


1.3679E-03 

1.3679E 

04 

4.1036E 

04 

I 

to 

288. 


1.2566E-03 

1.2566E 

04 

3.7699E 

04 

I 


302. 


1.154 8E— 03 

1.1548E 

04 

3.4643E 

04 

I 


316. 


1.0614E-C3 

1.0614E 

04 

3.1843E 

04 

I 


331. 


9.7589E-G4 

9.7589E 

03 

2.9277E 

04 

I 


347. 


3.9749E-C4 

8.9749E 

03 

2.6925E 

04 

I 


363. 


8.2561t-04 

8.2561E 

03 

2.4768E 

04 

I 


360. 


7.5968E— 04 

7.5968E 

03 

2.2791E 

04 

I 


398. 


6.9921E-G4 

6.9921E 

03 

2-0976E 

04 

I 


417. 


6.4372E-C4 

6.4372E 

C3 

1.9312E 

04 

I 


437. 


5.9280E— 04 

5.9280E 

03 

1.7784E 

04 

1 


457. 


5.46C4E-C4 

5.46G4E 

03 

1.6381E 

04 

I 


479. 


5.0310E-04 

5.0310E 

03 

1.5093E 

04 

I 


5C1. 


4* 6366H—04 

4.6366E 

03 

1.3910E 

04 

i 


525. 


4.2742E-C4 

4.2742E 

03 

1.2B23E 

04 

i 


550. 


3.94 1 1 E— 04 

3.9411E 

03 

U1823E 

04 

i 


575. 


3.6349E-C4 

3.6349E 

03 

1.0905E 

04 

i 


603. 


3.3533E— C4 

3.3533E 

03 

1.O060E 

04 

i 


631. 


3.C944E-C4 

3.0944E 

03 

9.28 3 IE 

oi 

i 




2.R561E-G4 

2.8561E 

03 

8.5683E 

03 

i 


692. 


2 . 6369E-04 

2.6369E 

03 

7.9106E 

03 

i 


724. 


2.4351E-C4 

2. 43 51 E 

03 

7.3052E 

03 

i 


759. 


2.2493E-C4 

2.2493E 

03 

6.7479E 

03 

i 


794. 


2.0783E-C4 

2.0783E 

03 

6.2348E 

03 

i 


832. 


1.9207E-04 

1.9207E 

03 

5.7622E 

03 

l 


671. 


1.7756E-C4 

1.7756E 

C3 

5.3269E 

03 

i 


912. 


1 .6420E-C4 

1.6420E 

03 

4.9260E 

03 

i 



Figure 
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955. 

1.5188E-C4 

1.5188E 

03 

4.5565E 

03 

I 

1000. 

1.4C54E-C4 

1.4054E 

03 

4.2161E 

03 

1 

1C47. 

1.3008E-04 

1.30C8E 

03 

3.9024E 

03 

I 

1C96. 

1.2C44E-C4 

1.2044E 

03 

3.6133E 

03 

I 

1148. 

1. 1 156E-C4 

1. 1156E 

03 

3.3469E 

03 

' I 

1202. 

1 .0338E-04 

1.0338E 

03 

3.1013E 

03 

I 

1259. 

9 • 5832E-05 

9.5832E 

02 

2.8750E 

03 

1 

1318. 

8.8878E-05 

8. 88 78 E 

02 

2.6664E 

03 

I 

1380. 

8.2469E-C5 

8.2469E 

02 

2.4741E 

03 

I 

1995. 

7.6562E-05 

7.6562E 

02 

2.2969E 

03 

I 

1514. 

7.111 9E— 05 

7.1119E 

02 

2.1336E 

03 

I 

1585. 

6.6102E-C5 

6.6102E 

02 

1.9831E 

03 

I 

1660. 

6 . 148CE-05 

6.1480E 

02 

1.8444E 

03 

I 

1738. 

5.7221E— 05 

5.7221E 

02 

1.7166E 

03 

I 

1820. 

5 . 3298E-05 

5.3298E 

02 

1.5990E 

03 

I 

1905. 

4.9685E— 05 

4.9685E 

02 

1.4906E 

03 

I 

1995. 

4.6348E-05 

4.6348E 

02 

1.3904E 

03 

11 

2C89. 

4.3253E-G5 

4.3253E 

02 

1.2976E 

03 

II 

2188. 

4.0384E-C5 

4.0384E 

02 

1.2115E 

03 

II 

2291. 

3.7723E-C5 

3.7723E 

02 

1.1317E 

03 

II 

2399. 

3. 5256E-C5 

3.5256E 

02 

1.0577E 

03 

II 

2512. 

3.2968E— 05 

3.2968E 

02 

9.8904E 

02 

II 

2630. 

3.0847E— 05 

3.0847E 

02 

9.2542E 

02 

II 

2754. 

2 . 8881 E— 05 

2.8881E 

02 

8.6642E 

02 

11 

2884. 

2.7057E-C5 

2.7057E 

02 

8.1171E 

02 

II 

3020. 

2.5365E— 05 

2.5365E 

02 

7.6096E 

02 

II 

3162. 

2 .3 796E-C5 

2.3796E 

02 

7. 1387E 

02 

II 

3311. 

2.2339E— 05 

2.2339E 

02 

6.7017E 

02 

II 

3467. 

2.0987E-C5 

2.0987E 

02 

6.2960E 

02 

II 

3631. 

1.9730E-05 

1.973CE 

02 

5.9191E 

02 

II 

3802. 

1 . 8563E-05 

1.8563E 

02 

5.5690E 

02 

11 

3981. 

1 . 7478E-C5 

1.7478E 

02 

5.2434E 

02 

II 

4169. 

1 . 6469E— 05 

1.6469E 

02 

4.9406E 

02 

II 

4365. 

1.5529E-05 

1.5529E 

02 

4.6587E 

02 

II 

4571. 

1.4654E-05 

1.4654E 

02 

4.3962E 

02 

II 

4786. 

1.3838E-05 

1.3838E 

02 

4.1514E 

02 

II 

5012. 

1.3077E-C5 

1.30 77E 

02 

3.9231E 

02 

II 

5248. 

1.2366E-05 

1.2366E 

02 

3.7G99E 

02 

II 

5495. 

1 . 17G2E— 05 

1.1702E 

02 

3.5107E 

02 

II 

5754. 

1.1081 E— 05 

1.108 IE 

02 

3.3244E 

02 

II 

6026. 

1.C500E-05 

1.05C0E 

02 

3.1499E 

02 

II 

6310. 

9.9550E-06 

9.9550E 

01 

2.9865E 

02 

II 

6607. 

9.4439E-06 

9.4439E 

01 

2.8332E 

02 

II 

6918. 

8.9641E— 06 

8.9641E 

01 

2.6892E 

02 

II 

7244. 

8.51326-06 

8.5132E 

01 

2.5540E 

02 

II 

7586. 

8.0890E— 06 

8.0890E 

01 

2.4267E 

02 

II 

7943. 

7.6895E— 06 

7.6895E 

01 

2.3069E 

02 

II 

8318. 

7.3130E-06 

7.3130E 

01 

2.1939E 

02 

II 

8710. 

6.9578E-06 

6.9578E 

01 

2.0873E 

02 

II 

9120. 

6.6223E-C6 

6.6223E 

01 

1.9867E 

02 

II 

9550. 

6.3053E-C6 

6.3053E 

01 

1.8916E 

02 

II 

icebo. 

6.C054E-C6 

6.C054E 

01 

1.8016E 

02 

II 


Figure 9 - Sample Output Listing-Program 5875 (Cont'd.) 
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SSLTUP LB4 UlSKrPLOTOliSAVE 

SEXECUTL It3 JOB 

IbJOB 

IBFTC MAiiN 

MAIN PROGRAM 3875 - LOCAL BUCKLING OF ISOTROPIC PANELS 

100 COMMON /INPUT/ AB(30) »BASlS(2> #BR(20) »BT(20) »CAPK(30) »CKNM(30> » 

1 LI t E2 1 F MU # IBUG # Ki, K2» K3# NAB# NBR # NOP# NRTlN, NBT# 

2 R T (200) # TITLL<12J# XLAB1(12>» YLAB1(12>» YLAB2<12># YLAB3(12)# 

3 YdSIG# YTSIG# YBBK# YTBR# YbSIGE# YTSIGE 

COMMON /COMPUT/ C» I tj # N# NRT # NRTC * OMIMU2# SIGRE# SIGPL# 
i SIOCRL ( 200 » 20 ) # S 1GCR 1 ( 200 #20 ) » SXGCR2 ( 2QQ # 20 ) 

150 SKIP = 0. 

200 WRITE ( 6 » 201 ) 

201 FORMAT ( 1 H 1 /// 49X » 34HLOC AL BuCKLINo OF ISOTROPIC PANELS ////} 

250 CALL uATAIN (SKIP) 

300 OM1MU2 = l.-FMU**2 
325 IT (Kl.Nt.l) GO TO 3b0 
026 NRT - NR TIN 
327 NRTC = 0 
340 GO TO 400 

350 CAi_l. PTC ( RT # NRT IN # NKT » NRTC ) 

CALL SETmIV ( 286 # 0 » 24 # 24 ) 

CALL SMXYV (1#1) 

400 GO iu ( b50» 550# 55U# 450# 450# 450» 450 )» NOP 
450 CALL CALC(l) 

500 IF ( NOP .LG. 4 ) GO TO 700 

501 IF l NOP .LQ. 5) GO TO 650 
550 CALL. CALC ( 2 ) 

600 IF (NUP.LG.2) GO TO 700 

601 IF ( NOP .EG. 6) GO TO 700 
6bU CAll CALL ( 3 ) 

700 CONllNUL 
720 GO TO 20u 
750 END 
SIBFTC bLKUAT 
llU BLOCK DA I A 

ill COMMON /INPUT/ AB(3U) »BASIS(2) »BR(20) #BT (20) »CAPK(30) #CKNM(30) # 

1 Li t L2 # FMiU » I BUG » Kl» K2# K3» NAB # NBR r NOP» NRT IN » NBT t 

2 K T 1 2 0 0 ) # T I ILL ( 12) » XLAB1U2)# YLAB1 ( 12 ) » YLAB2(12)# YLAB3(12)» 

3 YUiIG» YTSIG# YB8 k» YTBR# YbSlGE# YTSIGE 

120 DATA TITLE(l)# YLAblllJr YLaB 2(1), YLAB3(1)» XLABl(l) 

1 

2 / 36HL0CAL BUCKLING OF ISOTROPIC PANELS » 

3 36n BUCKLING STRESS# PSI FOR ALUMINUM # 

H 36HBUCKLING STRESS / ELASTIC MODULUS , 

5 36H BUCKLING STRESS# PSI FOR STEEL » 

6 36H RADIUS / THICKNESS / 

121 DA I A ( BASIS ( IB ) > 1B=1»2> /bH I #oH II /» 

* LI# L2 1 FMU / 1 • L7 » 3#L7# 0.30 / 

122 END 
SIBFTC DAT AIN 

SUBROUTINE DATAIN (SKIP) 

C SUBROUTINE TO READ IN AND PRINT OUT INPUT DATA. 
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TABLE VI - Fortran Listing-Program 3875 (Coni 1 d. ) 

C " ------ ~ 

50 COMMON /INPUT/ AB(30) # BASIS (2) #BR(20) »BT(20) »CAPK (30) ► CKNM (30) # 

1 Li# E2 » FMU# I BUG# Kl, K2, K3» NAB# NBR# NOp# NRTIN# NBT » 

2 RT ( 200 ) » TITLE ( 12) # XLAB1 ( 12) # YLAB1 ( 12) » YLAB2(12)» YLAB3(12)» 

3 TBSIG# YTSIG# YBBR# YTBR# YBSIGE, YTSIGE 

51 COMMON /COMPUT/ C# IB# N# NRT » NRTC # 0MIMU2# SIGRE# SIGPE# 

1 SIGCRE (200 # 20 ) » SIGCR1 (200 #20 ) # SIGCR2 ( 200 ► 20 ) 

75 IF (SKIP.NE.O.) GO TO 299 

100 READ (5»i0l) NAB# NBR# NBT# NRTIN# NOP# IBUG, Kl# K2» K3 

101 FORMAT (915) 

200 WRITE (6»201 ) NAB# NBR# NBT# NRTIN# NOP# IBUG# Kl# K2# K3 

201 FORMAT (4X»5HNAB =I5#4X,5HNBR =I5»4X#5HNBT = 15# 4X » 7HNRTIN =I5#4X» 
* 5HN0P =I5#4X»6HIBUG =I5»4X»5H Kl =I5»4X»5H K2 =I5»4X,5H K3 =15) 

250 IF (K2. EG. 0) GO TO 300 

251 SKIP = K2 - 1 

252 GO TO 300 

299 SKIP = SKIP - 1. 

300 READ ( 5 # 301 ) ( AB ( K ) » CKNM ( K ) » CAPK (K ) » K=1 # NAB ) 

301 FORMAT ( E10 . 5# A1 0 #E10 . 5 ) 

400 WRITt ( 6 # 401 ) ( AB ( K )» CKNM ( K) # CaPK ( K )» K = 1 # NAB) 

401 FORMAT ( /// ( 33X ► 20HPANEL ASPECT RATIO =F10 . 3# 16X » A10 » F10 . 3) ) 

500 READ ( 5# 501 ) (BRl J) # J=1#NBR) 

501 FORMAT (6E10.5) 

600 WRITE (6#601) (BR( J) , J=1 ,NBR) 

601 FORMAT ( /// 56X#2GHPANEL WIDTH / RADIUS //(5F23.6)) 

700 READ ( 5# 701 ) (RT(I) #1 = 1, NRTIN) 

701 FORMAT (6E10.5) 

800 WRITE ( 6 » 801 ) (RT ( I ) » 1=1 #NRTIN) 

801 FORMAT l /// 57X # 18HRADIUS / THICKNESS // (5F23.2) ) 

650 IF (Kl.EG.l) GO TO 1050 

900 READ ( 5# 901 ) YBSIG # YTSIG # YBBR » YTBR # YBSIGE# YTSIGE 

901 FORMAT ( 6E10.5 ) 

1000 WRITE (6# 1001 ) YBSIG# YTSIG# YbBR# YTBR# YBSIGE# YTSIGE 

1001 FORMAT ( // 4X» 6HY6S1G=1PE12.2» 3X» 6HYTSIG=1PE12.2» 

1 3X » 5hYBBR=lPE12.2# 3X# 5HYTBR=1PE12.2» 

2 3X# 7HYBSIGE=1PE12.2# 3X» 7HYTSIGE=1PE12 .2 ) 

1050 IF (NBT. EG. 0) GO TO 2000 

1100 READ ( 5# 1101 ) (BT(I) »I=1»NBT) 

1101 FORMAT (t.ElQ.5) 

1200 WRITE (6# 1201 ) < BT ( I ) # 1=1 » NBT) 

1201 FORMAT (/// 54X»23HPANEL WIDTH / THICKNESS // (6F23.5) ) 

2000 RETURN’ 

3000 END 
SlBFTC R1C 

SUBROUTINE RTC < RT » NRT IN#NRT #NRTC ) 

TO GENERATE RT(I)S FOR I FROM Nl TO N2» 

EVENLY SPACED RELATIVE TO LOGIC SCALE. 

DIMENSION RT ( 200 ) 

100 Nl = NRTIN + 1 
150 NRT = NRTIN + 101 
200 RT(Ni) = 100. 

250 RT (NRT ) =10000. 

300 DLOG = .02 
350 ONE =2.0 
400 NRTC = 101 
H50 00 650 1 = 1 » 99 
500 Nl = Nl + 1 

550 RTLOG = ONE + FLOAT ( I ) *DLOG 
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TABLE VI - Fortran Listing-Prograa 3875 (Confd. ) 

600 RT(N1) = 10 . **RTL0G 
650 CONTINUE 
700 RETURN 
750 END 
SlBFTC CALC 

SUBROUTINE CALC(NN) 

SUBROUTINE TO CALCULATE SIGMA CRITICAL AND SIGMA CRITICAL / E . 
SUBROUTINE NOVA IS CALLED TO MAKE PLOTS. 

10 COMMON /INPUT/ A6 < 30 ) , BAS IS (2 ) , BR ( 20 ) , BT < 20 ) , CAPK ( 30 ) , CKNM ( 30 ) , 

1 Ei, £2, FMU , IBUG, Kl, K2, K3, NAB, NBR » NOP, NRTIN, NBT, 

2 RT 1200)# T ITLE ( 12 ) , XLAB1112), YLAB1(12)» YLAB2C12), YLAB3U2), 

3 YbSIG, YTSIG , YBBR » YTBR * YBSlGE , YTSIGE 

11 COMMON /COMPOT/ C, IB. Nr NRT » NRTC. 0M1MU2. SIGRE. SIGPE. 

1 SIGCREC200.20) . S IGCR1 ( 200 . 20 ) . S1GCR2 ( 200 . 20 ) 

N — N N 

50 DO 1090 K=lrNA8 

70 DO 1030 u=lr'NBR 

71 IF (bK< J) .NE.0.0) 60 TO 90 

72 wRiTt (6.73) N. AB ( K } » CKNM(K) . CaPK(K). BR(J) 

73 FORMAT ( 1H1 , 62X , 7HOPJ ION ,11 /// 14X. 19HPANEL ASPECT RATI0=F8.2, 

1 14X,A10,F8.3»l5X, 19HPANEL WiDTH/RADIUS=FlO .4 /// 

2 11X.21HPANEL WIDTH/THICKNESS 9X, 17HBUCKLING STRESS/E 13X, 

3 19HBUCKLING STRESS, PSI 11X» 19HBUCKLING STRESS, PSI/ 

4 73X.14HL = 10,000,000 16X.14HE = 30,000,000 /// ) 

74 DO 64 1=1, NBT 

75 IF (bT(I) .EQ.0.0) GO TO 82 

76 SIOCKE(I.J) = CAPK(K) / OMIMU2 * ( 1 . /BT ( I ) ) **2 

77 SIGCR1 ( I , J) = SIGCRE(I.J) * El 

78 SI6CR2 ( I , J) = S1GCRE ( I , J) * E2 

79 WRITE (6,80) BT(I), SlGCRE<IrJ)r SIGCRl(I.J), SIGCR2(I,J) 

80 FORMAT (1H F25.4, 3(1PE30.4>) 

81 GO TO 84 

82 tfRllt ( 6 r 83 ) BT ( I ) 

83 FORMAT ( 1H F25.4, 1UX, 31H CALCULATIONS OMITTED, B/T = 0.) 

64 CONTINUE 

65 GO Tu 1030 

90 WRITE (6,91) N, AB (K) , CKNM(K), CAPK(K), BR(J) 

91 FORMAT ( 1 Pi 1 , 62 X , 7H0PT I ON ,11 /// 14X , 19HPANEL ASPECT RATI0=F8.2, 

1 I4X,AI0,F8.3»15X,19HPANEL WIDTH/RADIUS=F10 .4 /// 

2 7X»I6HRADIUS/THICKNESS,11X,17HBUCKLING STRESS/E » 12X , 

3 19HBUCKLING STRESS, PSI , 11X, 19HBUCKLING STRESS, PSI .6X.5HBASIS / 

4 t>6X,l4h£ = 10,000,000 16X,l4HE = 30,000,000 // ) 

92 NOP = NOP 
100 00 1020 1=1, NRT 

120 C=U.oG5-(0.546*(l.-EXP{-l./lb.*SGRT(RT(I) ) ) ) ) 

130 GO TO (150,210,350) ,N 
150 S 1 ORE = C* ( 1 «/RT ( 1 ) ) 

190 GO 10 490 

210 IF (AB(K)*BR( J) .GE.1.0) Go TO 290 

230 SIGRE = 11.28*RT(I)**(-l.b) + 0 . 109* ( RT ( I ) *AB ( K ) *BR ( J) ) ** (-1 . 3 ) 

270 GO TO 490 

290 SIGRE = 11.26*RT(I)**(-1.6) + 0 . 109* ( RT ( I ) *AB ( K ) *BR ( J ) ) ** ( -1 . 3 ) 

* - 1.416 * RT ( 1 ) ** (“1 .6 ) * ALOG ( AB ( K ) *BR ( J ) ) 

330 GO TO 490 

350 IF (AB(K)*BR(J) .GE.1.0) GO TO 430 

370 SIGRE = o.0l*RT(I)**(-1.6) + 0 .076* (RT ( I ) *AB(K >*BR ( J) ) ** (-1 .3) 

410 GO TO 490 
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430 SlGKt = 7.52 * RT ( I ) ** < -1 .6 > + 0.0?2 * l RT < I ) *AB (K > *BR < J ) ) ** ( -1 . 3 ) 
* - 1.418 *RT ( I ) **(-i«6) * ALOg ( AB ( K ) *BR (U> ) 

490 SIGPE = CAPK(K) / 0M1MU2 * ( 1 . / ( RT ( I ) **2 * BR(J)**2>) 
o30 IF (SIGRE .LE. 2.*S1GPE) GO TO 730 
650 IB = 2 

670 SIGCRE ( I » J) = SIGRE 
710 GO TO 770 
730 IB = 1 

750 5IGCRE ( I » J) = SIGPE + SIGRE**2 / (4. * SIGPE) 

770 CONTINUE 

760 SlGCkl { I * JJ = SIGCREI I ,<J) * El 
790 SIGCR2U,J> = SIGCRE 1 1 » J) * E2 
610 IF l IBUG.NE* 0 > GO TO 1010 
630 I.F (I.LE.NRTIN) GO TO 1010 
650 GO TO 1020 

1010 WR II L (6» 1011 ) RT(I> »SIGCRE(1,J> ,SIGCR1(I,J> ,SIGCR2(I,J) ,BAS1SIIB) 

1011 FORMAT (1H r F17 . 0 , 3 ( 1PE30 . 4 ) , Ai6 ) 

1020 CONTINUE 

1030 CONTINUE 

1050 IF (Kl.EQ.l) GO TO 1090 
1070 CALL NOVA(K) 

1090 CONTINUE 
2010 RETURN 
2030 END 
$I6FTC NOVA 

SUBROUTINE NOVA (K) 

COMMON /INPUT/ AB(3G) ,BASIS(2) ,BR (20 ) ,BT ( 20 ) , CAPK (30 ) ,CKNM(30> , 

1 El, E2, F MU . IBUG, Kl. K2, K3, NAB , NBR, NOP, NRTIN, NBT » 

2 RT(20U), T ITLE ( 12 ) » XLAB1 ( 12 ) , YLaBI ( 12) , YLAB2U2), YLAB3Q2), 

3 YuSIG, YTSIG, YBBR. YTBR, YbSIGE, YTSIGE 
COMMON /COMPUT/ C, I8r N, NRT , NKTC , 0MIMU2, SIGRE, SIGPE, 

1 SiGCRE(200,20) . S1GCR1 ( 200 , 20 ) , SIGCR2 ( 200 ,20 ) 

TO PLOT THREE FRAMES WITH NBR CURVES EACH FRAME 
LOO-LOG SCALE FOR ALL PLOTS 

PLOT 1 SIGCR1 v S R/T ( E=10 , 000 , 000 ) 

PLOT 2 SIGCR2 VS R/T ( E=30 , 000 , 000 ) 

PLOT 3 SIGCRE V S R/T 

CALL SETMIV ( 288 ,0,24,60) 


CALL SUBPLT (SIGCR1, YLABl , YBSIG, YTSIG, K ) 

CALL SUBPLT (SIGCR2, YLAB3* YBSIG, YTSIG, K ) 

CALL SUBPLT (SIGCRE, YLAB2 » YBSIGE, YTSIGE, K ) 

RETURN 

END 

SIBFTC SUBPLT 

SUBROUTINE SUBPLT (F, YNAME , YB, YT ,K) 

COMMON /INPUT/ AB<30) , BASIS <2> , BR <20 ) ,BT(20) , CAPK (30 ) ,CKNM(30) , 

1 El, E2 , FMU, IBUG, Kl, K2, K3, NAB, NBR, NOP, NRTIN, NBT, 

2 RT (200) , TITLE ( 12 ) , XLAB1 ( 12) , YLABl ( 12 ) , YLAB2 ( 12) , YLAB3(12), 

3 YBSIG, YTSIG, YBBR, YTBR, YbSlGE, YTSIGE 

COMMON /COMPUT/ C, IB, N, NRT, NRTC, 0MIMU2, SIGRE, SIGPE, 

1 SIGCRE (200,20) , SIGCR1 ( 200 ,20 ) , SIGCR2 (200 , 20 ) 

DIMENSION OPT (2,3) 
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DIMENSION F ( 200 » 20) r YNAME ( 12) 

DATA (OPT ( 1 r I ) r 1 = 1 r 3) / 12H OPTION 1 r 12H OPTION 2 
* 12h OPTION 3 / 


AA= Ad ( K ) 
CC=CaPK £ K ) 

102 XR = 1.E4 

103 XL =100. 


i 


N=N 

CALL GRIDlV(4rXLrXRr YBrYTr 1.# 1., Or Or Or Or 5r 6) 

C PRINT TOP TITLEr VERT r HORIZ 

CALL APRNTV (Or -16r 36r YNAM'Er 200r 900) 

CALL PRINTV ( 36 » XLAblr &00r 12) 

CALL PRINTV (-6r 5HA/b.= r300r 980 .) 

CALL PRINT V ( -6 r SH K = rfcOOr yQO ) 

CALL PRINTV ( 12 r OPT ( 1 rN) r900r 980 ) 

CALL LABLV ( AAr3o0r980r 6» 1»3> 

CA.-L L ABC V (CCr 660 r 980 r 6 » 1 » 2) 

CALL kITc2v (380rl0lLri023r90rir 36 r -1 r 36HL0CAL BUCKLING OF ISOT 
*ROP IC PANELS r NLAST ) 

DO i5U J=1 r NbR 
Nb=NRT IN+ 1 

IF IbK(J)) I50rl50rlul 
101 DO luU I=NB r NRT 

IF (F ( I fo) .LE. YT) oO TO 10b 
100 CONTINUE 
GO TO ISO 
10b Nb=i+1 

DO i Ob 1=NB r NRT 
1Y = NRT+N8 -I 

IF (F(1Y»J) .GE. YB) GO TO 110 
106 CONTINUE 
GO Tu 150 
110 NL=1 Y 

kT i = RT(Nb-l) 

FF = F (Nb-lrJ) 

1X1 = NXV(RTT) 
i Y 1 = NY V ( FF) 

DO 140 l=NBrNE 
KIT- KT ( I ) 

FF= F(lrJ) 

1X2 = NXV(RTT) 

IY2 = NYV ( FF) 

DO 130 NF= lrb 

CALL LINEV (IXlr IYlr IX2» IY2) 

130 CONTINUE 
IX1=1X2 
1Y1=1Y2 
140 CONTINUE 
150 CONTINUE 
200 RETURN 

end 

12 END OF FILE OR DATA CARD DATA FOLLOWS 

4 7 0 2 4 1 

.4 KS=KC= 7.0 

•6 KS= 4.0 

.6 KC= 5.7 

.8 KS= 3.29 
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.01 .0 2 .05 .1 

i # 

100 . 1000 . 

100 . 100000 . 


.2 .5 

.ooooi .0* 
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APPENDIX A 

SUPPLEMENTARY OPTIONS 


A. 1 OPTION 2 

A. 1.1 GENERAL 

All of the curves presented in Figure 11 are based upon substitution 
of the OPTION 2 full cylinder equations [see equations (2-9) and (2-10)] 
into the Schapitz criterion [2]. These curves were obtained by using 
digital computer program 3875 and an automatic plotter. The program is 
described in SECTION 5 and can be used to obtain additional plots and/or 
single-point solutions, as desired. Since the automatic plotting machine 
does not provide any capability for the print-out of lower case letters, 
the ratio (a/b) appears on all the given buckling curves as (A/B). In 
view of the fact that the OPTION 2 full cylinder equations constitute a 
mean-expected criterion, the curves given here may be of value in the 
comparison of predictions versus test data. In addition, when used in con 
junction with the buckling curves based on OPTIONS 1 or 3, these plots 
give an indication of the dispersion which might be expected among a large 
array of test data. 

A. 1.2 BUCKLING CURVES 

Table VII lists the families provided here. All of these curves are 
based on elastic behavior and a Poisson's ratio of 0.30. 
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A. 2 OPTION 3 

A. 2.1 GENERAL 

All of the curves presented in Figure 12 are based upon substitution 
of the OPTION 3 full cylinder equations [see equations (2-ll) and (2-12)1 
into the Schapitz criterion [2], These curves were obtained by using digital 
computer program 3875 and an automatic plotter. The program is described 
in SECTION 5 and can be used to obtain additional plots and/or single-point 
solutions, as desired. Since the automatic plotting machine does not provide 
any capability for the print-out of lower case letters, the ratio (a/b) 
appears on all the given buckling curves as (A/B). It is recommended that 
these plots be used whenever (a/R) > 1 since, in this range, the curves of 
SECTION 4 can give unconservative design values. 

A . 2 . 2 BUCKLING CURVES 

Table VIII lists the families provided here. All of these curves are 
based on elastic behavior and a Poisson's ratio of 0.30. 


A-18 


GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 



TABLE 

VIII - Table 

of Contents for the 




Supplementary Buckling 
Curves Based on OPTION 3 


Figure 

(a/b) 

K 

Page 

12 

(a) 

.4 

7.00 

A-20 

12 

(b) 

.6 

4.00 

A-21 

12 

(c) 

.6 

5.70 

4 A-22 

12 

(d) 

.8 

3.29 

A-23 

12 

(e) 

.8 

5.70 

A-24 

12 

(f) 

1.0 

3.29 

A- 25 

12 

(g) 

1.0 

5.70 

A-26 

12 

(h) 

2.0 

3.29 

A-27 

12 

(i) 

2.0 

5.70 

A- 28 

12 

(j) 

4.0 

3.29 

A^29 

12 

(k) 

4.0 

5.70 

A-30 

12 

(1) 

10.0 

3.29 

A-31 

12 

(m) 

10.0 

5.70 

A-32 

12 

(n) 

100.0 

3.29 

A-33 

12 

(o) 

100.0 

5.70 

A-34 


A- 19 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 


B 

U 

c 

K 

L 

I 

N 

6 

S 

T 

R 

C 

3 

S 

/ 

e 

L 

A 

S 

T 

t 

C 

M 

O 

0 

0 

L 

U 

* 


BUCKLINS OF ISOTROPIC PANELS 

A/t.a 0.40 K a 7.000 0RT10R 5 



Ar-20 


GENERAL DYNAMICS CONVAIR DIVISION 


GDC-DDG-67-006 




* A C 9 


GDC-DDG-6 7-006 



RADIUS / THICKNESS 


Figure 12(c) 


A.-22 


GENERAL DYNAMICS CONVAIR DIVISION 


U1 «J << *> »~ 


GDC-DDG-67-006 



>000 

RADIUS / THICKNESS 


Figure 12(d) 
Av-23 


GENERAL DYNAMICS CONVAIR DIVISION 


-4 o* > r n 


GDC- DDG-6 7-006 


BUCKLING OF ISOTROPIC PANELS 

A/ta« 04*0 K « 5 4 700 OPTION 3 



A^24 

GENERAL DYNAMICS CONVAIR DIVISION 


GDC-DDG-67-006 



1000 

RAOIUf / THlCKNCtS 


Figure 12(f) 


A-25 


GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 


6 

U 

c 

K 

L 

I 

N 

6 

S 

T 

R 

E 

» 

S 

/ 

E 

L 

A 

S 

T 

l 

C 

M 

O 

0 

U 

L 

U 

• 


BUCKLING OF ISOTROPIC PANELS 

A'!** 1*00 K s 5.700 OPTION S 



A— 26 


GENERAL DYNAMICS CONVAIR DIVISION 


h e» ► i” Ft v •> v> m » -t w 


GDC- DDG- 6 7-006 


BUCKLING OF ISOTROPIC PANELS 

*/».* 2.00 K = 3 >290 OfTlOM 3 



A-27 

GENERAL DYNAMICS CONVAIR DIVISION 



h «i > r n n in ti n iv *1 (• 


GDC-DDG-6 7-006 


B 

u 

c 

K 

L 

t 

N 

6 


1 

C 

M 

O 

D 

U 

L 

U 

B 


BUCKLING OF ISOTROPIC PANELS 


*/B.» *.00 K a 5.700 OPTION 5 

. 0*000 



A— 28 

GENERAL DYNAMICS CONVAIR DIVISION 


h w » r n v w w pi » -» v *> X 


GDC- DDG- 6 7-006 


B 

U 

C 

K 

L 

I 


I 

C 

M 

O 

0 

U 

L 

U 

s 


*/».■ 4.00 

.otooo 


BUCKLING OF ISOTROPIC PANELS 

K * 3.290 OPTION 3 


Figure 12 (,j) 
As-29 



GENERAL DYNAMICS CONVAIR DIVISION 


GDC-DDG-67-006 


BUCKLING OF ISOTROPIC PANELS 

4.00 K « 5.700 OPTION 3 



A-30 


GENERAL DYNAMICS CONVAIR DIVISION 



GDC-i)iX5- 6 7-006 


B 

U 

c 

K 

L 

t 

N 

6 

3 

T 

R 

C 

s 

3 


E 

L 

A 

3 

T 

I 

C 

M 

O 

D 

U 

L 

U 

3 


BUCKLING OF ISOTROPIC PANELS 

A/*«» 10.00 K s 3.290 OPTION 3 



A-31 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC-DDG-6 7-006 


BUCKLING OF ISOTROPIC PANELS 


A/t.a 10.00 K s 1.700 OPTION 3 



A*-32 

GENERAL DYNAMICS CONVAIR DIVISION 


< ‘ 


GDC- DDG-6 7-006 





GDC- DDG-6 7-006 


B 

U 

C 

K 

L 

1 

N 

0 

S 

T 

R 

e 

s 

s 

/ 

e 

L 

A 

s 

T 

1 

C 

M 

O 

D 

U 

L 

U 

» 


BUCKLING OF ISOTROPIC PANELS 

*/».» 100.00 K s 9. TOO OPTION S 



A— 34 

GENERAL DYNAMICS CONVAIR DIVISION 


